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ABSTRACT 


The  one-on-one  stochastic  duel  is  extended  to  tun  general  two-on-one  duel 
for  the  first  time.  The  state  equations,  win  probabilities,  means  value  and 
variance  functions  are  derived.  The  case  where  one  side  has  Erlang  (2)  fir¬ 
ing  times  and  the  other  is  negative  exponential  is  compared  to  the  corres¬ 
ponding  "Stochastic  Lanchester"  and  Lanchester  models  to  demonstrate  their 
non-equivalence. 


*  Th i s  is  an  expanded  version  of  a  report  by  the  same  title  published  as 
ISE  TR  83-1  by  the  Department  of  Industrial  and  Systems  Engineering,  Univer¬ 
sity  of  Southern  California,  Los  Angeles,  CA  70089 


1. 


INTRODUCTION 


In  this  paper  we  extend  the  theory  of  stochastic  duels  to  the 
two~on-one  case.  Previously  (see  Ancker  (1)  )  much  work  has  been  done 
on  the  marksman  problem  (one  versus  a  passive  target)  and  the  one-on- 
one  duel  but,  to  our  knowledge,  this  is  the  first  general  extension 
which  does  not  necessarily  assume  all  interfiring  times  have  the 
negative  exponential  distribution  (ned)  (Clark  (4])  or  are  constant 
(Ancker  (2],  Anderson  [3]). 

In  the  following  sections  the  problem  is  rigorously  formulated, 
and  the  state  equations  are  derived  and  solved.  Then  we  first  assume 
all  ned  interfiring  times  as  one  example  (this  checks  with  previously 
derived  "Stochastic  Lanchester"  equations)  and  secondly  show  an 
example  with  Erlang  (2)  interfiring  times  on  the  A  side  and  ned 
interfiring  times  on  the  B  side. 

These  two  examples  are  then  compared  with  the  equivalent 
Lanchester  Square  Law  solution  to  show  that,  for  such  small  numbers 
involved,  both  the  Lanchester  and  the  "Stochastic  Lanchester” 
approximations  to  a  situation  in  which  interfiring  times  are  not  all 


□  □ 


the  other  side  Is  destroyed.  There  are  two  contestants  on  the  A  side, 
who  fire  continuously  and  independently  of  each  other  and  who  have 
identical  random  interfiring  times,  XA ,  (with  probability  density 


function  (pdf)  f x  (x))  and  identical  kill  probabilities,  pA,  on  each 
round  fired.  The  B  side  has  only  one  contestant  whose  random 
interfiring  time  is  Xg  with  pdf  fx  (x)  and  kill  probability  pR.  All 
interfiring  times  are  independent.  Both  sides  start  at  time  zero  and 
each  contestant  fires  his  first  round  (if  still  alive)  at  one  random 
interfiring  time  later.  There  are  no  time  or  ammunition  limitations 
on  the  contest. 


3.  GENERAL  SOLUTION 


Our  solution  technique  has  three  principal  features.  One,  we 
shall  first  consider  each  contestant  separately  as  a  marksman  firing 
at  a  passive  target,  and  two,  we  shall  concentrate  only  on  firings 
which  are  kills  rather  than  on  every  firing  event.  This  is  not 
necessary  but  simplifies  the  mathematics.  This  procedure  entails 
knowing  the  solution  to  the  problem  of  the  marksman  firing  at  a 
passive  target  which  is  exhaustively  examined  in  reference  [1]  and 
which,  we  shall  consider  known.  Thus,  if  TA  is  the  time  to  a  kill  by 
one  of  the  A's,  we  shall  designate  the  distribution  function  (df)  and 
the  pdf  of  T^  as  FA(t)  and  fA(t)  respectively.  Similarly  for  B,  we 
have  Fg(t)  and  fg(t). 

Finally,  we  shall  superpose  the  three  marksmen's  firing  sequences 
and  use  the  backward  recurrence  time  technique  to  write  the  state 
probability  equations.  This  means  that  we  observe  the  superposed 


«>  i*cj  iinire  at  the  time  t  and,  if  we  define  Y  to  be  the  random  variable 
I  I  in  >:  i  ii  <•  <>  tin*  last  event  (kill)  t  h  e  n  the  first  or  d  e  r  p  r  <•  ha  h  i  1  i  t  y  that 
the  next  event  will  be  a  kill  by  A  and  will  occur  in  the  interval 
(t,  t  +  A)  is  given  by 


fA(y) 

r  (y)A  =  -  A  +  o(A) 

F'(y) 


for  A,  and  by 


fB(y) 

r  (y  )A  - - A  +  o(A  ) 

F^(y) 


(1) 


(2) 


c  c 

for  B.  FA(y)  an<^  Fjj(y)  are  t*le  complementary  distribution  functions 
(cdf's)  for  the  killing  times  of  A  and  B  respectively.  See  reference 
[1]  for  a  discussion  of  the  backward  recurrence  time  technique.  The 
r(y)'s  are  the  instantaneous  kill  rates  for  each  marksman.  We  note 
that 


fA(y) 
fa  (y ) 


(3) 


and  similarly  for  B.  We  note  that  if  there  have  been  no  events  at 
time  t,  then  y  =  t.  The  state  random  variables  are  defined  to  be. 


N^(t)  =  the  number  alive  on  side  A  at  time  t,  J 

>  (A) 

N^(t)  =  the  number  alive  on  side  B  at  time  t,  j 


and  the  state  proba  ,’lities  are  defined  as 


P13(t> 


P[i  are  alive  on  side  A,  and  j  are  alive 
on  side  B  at  time  t]  , 


and 

*’ij  (t,y  )dy 


P[i  are  alive  on  side  A,  j  are  alive  on 
side  B  at  time  t  and  it  has  been  between 
y  and  y  +  dy  time  units  since  the  last  event], 


i  -  0,1,2 
j  3  0,1. 


(5) 


Th»  relation  between  the  two  types  of  probabilities  in  (5)  is 
pij*  '  =  /q  pij  The  second  form  in  equation  (5)  is  necessary 

since  in  general,  at  any  arbitrary  time  t,  the  state  of  the  system  is 
defined  by  the  number  alive  on  each  sid.-  and  the  time  to  the  last 
event.  The  backward  recurrence  time  does  not  exist  for  the  beginning 
st'te  or  terminal  states.  In  fact,  for  our  small  system  it  is  only 
needed  for  state  N^(t)  *  1  and  Ng(t)  =  1,  since  at  most  two  events 
te.rainate  the  system.  In  more  complicated  systems,  more  backward 
recurrence  times  must  be  accounted  for,  but  we  only  need  one  since  at 
most  one  event  can  precede  termination. 

The  problem  will  be  to  derive  the  state  probabilities  and  from 
these  to  obtain  the  win  probabilities,  P [ A ]  and  P[B],  and  the  means 
and  variances  of  the  state  random  variables. 

The  state  probabilities  are  obtained  in  the  usual  fashion.  For 
example,  in  the  first  equation  in  (6)  below,  the  probability  of  being 


-A- 


In  state  NA(t)  =  2  and  t )  =  1  at  t  twe  t  +  A  is  equal  to  the 

probability  of  being  in  the  state  at  time  t  times  the  probability  of 
remaining  there  in  the  interval  (t,  t  +  A),  which  is  the  probability 
that  one  A  fails  to  kill,  (1  -  rA(t)A),  times  the  probability  the 
second  A  fails,  (1  -  r.(t)A),  times  the  probability  the  B  contestant 

A 

fails,  (1  -  r  (t)A).  All  these  probabilities  are  independent  by  the 
B 

statement  of  the  problem.  Hence, 


P21(t+A)  =  p21(t)  tl  -  rA(t)A]2  (1  -  r  g  ( t  )A  j  +  o  (A  ) 
and  s i mi lar ly , 

P2Q(t+A)  =  p2Q(t)  +  p21(t)  2r  A  ( t  )  A  (I  -  rfi(t)A]  +  o(A), 

Pn(  t+A  ,  y+A  )  =  pn(t,y)  [1  -  r  A  ( t )  A  )  [  1  -  rfi(y)A)  +  o(A), 

P1Q(t+A)  =  p1Q(t)  +  /o  P1 1  Ct  ’ y  5  1 1  "  ]rA(t)Ady  +  o(A), 

pQ1 (t+A  )  =  pQ1  (t  )  +  /o  P11  (t , y )rfi (y )Ady  +  o(A  ). 


(6) 


Rearranging  terms,  dividing  by  A  and  letting  A  +  0,  we  have  that 


dp2],(t) 

_ 

dp2Q(t) 

dt 


-  p  2  j  ( t  )  1 2  r  A  ( t  )  +  rg(t)3, 

2rA(t)  P21U), 


(7) 

(8) 


(IT  +  h  +  rA(°  +  rB(y)^  Pil(t>y)  =  °» 


11 


dpio(t) 

dt 


(t)  /q  pu(t,y)dy  =  rA(t)p,,(t), 


1 1 


dp01  >  r 

- =  / o  Pn  (t,y)rR(y)<fy, 


B 


(9) 

(10) 


dt 


(11  ) 


with  initial  conditions, 


P21(0)  -  1, 

P± j (0)  »  0  i  *  2,  j  *  1, 


(12) 


and  boundary  condition 


P  j  1  ( t  +  A,  0  +  A)  A  -  p21(t)rB(t)A  +  o  (  A)  , 

which,  on  dividing  by  A  and  letting  A  -*■  0,  gives 


Pn(t,0)  -  p21(t)rB(t). 


(13) 


Equations  (7)  and  (8)  are  easily  solved,  in  sequence,  using 
initial  condition  (12),  to  give 

P21(t)  "  1  FA(t)]2  FB(t)  (14) 

and 

p20(t)  “  2  /q  FAU)  FBU)  fAU)d5*  (15) 

The  general  solution  to  equation  (9)  is 


g(t-y)e  e 

Using  boundary  condition  (13), 


-  g(t-y) 
we  find 


Fd(y)  P'(t>. 


-6- 


(16) 


g(t-y)  =  F^(t-y )ffi(t-y) 
and  so  , 

Pu(t,y)  =  Fj(t)  F^(t-y)Fg(y)fg(t-y). 

Integrating  out  y  from  (t,y), 

P 1 1  Ct )  =  Fj(t)  /*F^(t-y)  FB(y)  fB(t-y)dy,  (17) 

which  satisfies  initial  conditions  (12).  Substituting  (16)  and  (17), 
into  equations  (10  and  (11),  we  have  at  once  that, 

p10(t)  =  /  ofA(UdC  /JFA(?_y)FB(y)fBU'y)dy>  (18) 

and 

P01(t)  =  ^oFACOdC  /oFAU_y)fB(y)fBU~y)dy  ’  (19) 

which  also  satisfy  initial  conditions  (12). 

Equations  (14),  (15),  (17),  (18)  and  (19)  are  the  equations  of 
state  for  the  model  and  sum  to  one  for  all  t  and  along  with  iMju.it  ion 
(16)  contain  essentially  all  the  information  on  the  solution.  Some 
other  expressions  of  interest  may  be  derived  as  follows. 

The  probabilities  that  A  wins,  P [ A 1  and  B  wins,  P [ B  j  are 

P[A)  =  P2Q  (»)  +  P1Q(~)  -  2  J*  F®(U  fAU)dC 

+  J"  fA(OdC  /*F*($-y)  FA(y)  fBU,y)dy  (20) 

-7- 


P  t  B  ]  =  p01(«)  =-  /“  F®(C)dC  l\  F  =  (C-y)  fB(y)  fB(C-y)dy, 
and  P[A]  +  P [ B ]  *  1,  as  it  should. 

The  means  and  variances  may  be  calculated  as  follows: 


E  [  NA  ( t )  ]  =  nA(t)  s  2^P2i^^  +  p20^t^  +  1^P11  +  PiQ^^’ 


E  [  Nfi  ( t )  1  =  nfi  ( t  )  »  l(p21(t)  +  Pn(t)  +  P01(t>), 


E I NA  <  t )  1  =  n2(t)  “  4(p21(t)  +  p2Q(t))  +  l(pn(t)  +  p1Q(t)), 


E[Ng(t)j  =  rig(t)  -  l(p21(t)  +  Pn(t)  +  pQ1(t))  ®  ng(t)  , 


VlNA(t)1  =  °A(t)  "  lnA(t)]2 
and 


V  t  N  B  ( t )  ]  =  nB(t)  -  [ng(t  )  ]  2  -  n  B  ( t )  -  UB(t)]\ 

The  above  are  the  general  results  for  the  stated  problem 


Examples  for  specific  interfiring  times  follow. 


4. 


K  XAMPLKS 


K  x  ;un  p  1  e  1 

Let  the  Interfiring  times  be 


Vx>  *  \ e 


,  f  v  (x)  =  —  e 


From  prior  work,  reference  [1] 


,  ,  ,  PA  PA 

f  .  (t  )  «  —  e 

A  “a 


...  PB  yB 
•  fBU)  '  ^  6 


PA  PB 

and  from  equations  (1)  and  (2)  r.  =  —  ,  r„  =  —  and  neither  is  a 

A  “a  b  pb 

function  of  time  to  the  last  event  (this  is  a  result  of  the  no  memory 


property  of  the  negative  exponential  random  variable). 

Using  f.(t),  fR(t)  and  their  cdf's  in  (14),  (15),  (17),  (18)  and 


(19)  we  have  easily  that, 


P21U> 


(rA  +  r  B^  P 


P  (  t  ) 
20 


'a  ,,  -<2rA+rB)C. 

-  [ 1-e  ]  , 


2r . +r „ 
A  B 


pu(t; 


rB  .  (rA+rB)t  ■(2rA+tB)t, 
—  (e  -  e  ]  , 

A 


pio(t) 


poi(t) 


rB  "'rA+I‘B)t,  rB  .  (2rA  +  rB)C  . 

vV*  ”  1  ’  * 1 1  ■ e  ’• 

rB _  ,,  -(rA+rB)t.  _ rl _ ,,  -(2rA+rB),; 

V'a+'b*  6  rA(2tA+rB) 


These  results  may  also  be  obtained  from  rhe  so-called  "Stochastic. 
Lanchester"  equations  which  are  a  stochastic  version  of  the  famous 
deterministic  Lanchester  equations.  In  particular  these  are  the  two- 
on-one  version  of  the  "Square  Law".  For  an  exposition  on  these 
matters  see,  for  example,  Clark  [A],  The  state  equat  ons 
corresponding  to  (7)  through  (12)  are  quite  similar  but  only  involve 
ordinary  differential  equations  and  are  somewhat  simpler.  We  also 
have  the  following  results  from  (20)  and  (21): 


P  [A] 


rA(2rA+3rB) _ 

(rA+rB)(2rA+rB)  ’ 


(29) 


P(B]  = 


(rA+rB)(2rA+rB)  ’ 


(30) 


and  from  (22)  and  (23), 


nA(t) 


rA(5rB+ArA) 


-(rA+rB),: 


(rA+rB)(2rA+rB)  rA(rA+rB) 


.  rB(rA‘rB)  “(2rA+rB}t 
rA(2rA+rB>  6 


(31) 


nB(t)  =  (rA+rB)(2rA+rB) 


“(rA+rB)t: 


(rA+rB} 


(2rA-rB)  -(2rA+rB)t 

<2VrB>  6 


(32) 


-10- 


and  f roa  (24)  and  (25), 


l1 ‘l. 

rA+rB)(2rA+rB) 


'(rA+rBU 


rB(3rA~rB 

rA(2rA+rB) 


-(2rA+rR)t 


(33) 


n  B ( C )  =  nB(t)  * 


(34) 


The  variances  may  now  be  calculated  using  equations  (26)  and  (27) 


Example  2 
Let 


f  Y  (X) 
*A 


f 


Again  from  reference  1 1 ] »  setting  qA  -  1“PA> 


In  this  case,  since  each  contestant  on  side  A  has  an  Erlang  (2) 
interfiring  time  (with  the  same  mean,  p^,  as  in  example  I),  r^< y)  is 

no  longer  a  constant  but  is  a  function  of  y,  -he  backward  recurrence 
time.  Proceeding  as  before,  after  much  integration  and  algebra  we 
have  the  following  results.  Using  the  notation 


»i  ’  h  (1  -  /i*  >•  a2  '  h  a  +  /q*  >•  rB  '  ST  • 

A  A  B 


then. 


rBt  a , t  -a_  t  _ 

.  .  e  ,  1  2  ,  Z 

P2l(t>  -  - - -j-  l«2*  -«!*  )  • 

(  a  2  -a  j  ) 


2aiao  i  -(2a. +i*  )t  . 

P20<t)  *  7«7^7  {a1(2a1+rB>  ll_e  1  +  vw  u“e 


-( 2a2+rg  )t 


(al'hx2) 

a  la2  (a  l+a2+rB) 


-(“2+“2+^B)t]  . 
ll  -  e  }  , 


90  99 

r_e  (a  „  -a ,  )  -a .  t  (a. -a,)  -a„t 


- (a , +a„ ) t 


P10(t) 


B 


.<02  '  °1>  .. 

_  [  l  -e  1 

(a  -a  )2  (al+rB) 


(a2"al}  ,,  _(a2  +  rB)t, 

(a 2  +b )~  (1'e  1 


(  2ot  1  +rB  ) 


-(2a  +rR)t 
[1-e  1  ®  1  - 


<2«2+rB) 


-  (  2a_+r  )  t 

u-«  2  B  i 


(«2  +  «2>,, 

+  U-e  ]} , 

(ai+a2+rB) 


(38) 


rB  ,  < a  2 _a  1  ^  _(al+rB^t: 

p_,(t)  - - r-  { - t - ~ — r  (1-e 

01  («2-ai)2  a  1  (a  I  "^rB  5 


1  _  r  '(n2+rB)t: 

a2(a2+rB) 


al(2a1+rB) 


-(2a1+rB)t  cxj 

[1-e  ]  - 


a2(2ot2+rB) 


-( 2a„+r  )t 
(1-e  2  B  ] 


(ol+a2+rB) 


-(a. +a0+r  )t 
(1-e  1  2  B  ]}, 


(39) 


PfAl  -  4ala2(al+q2)~2(al+a2)rB~(al+a2)rB  . 
(2a1+rB)(2a2+rB)(a1+a2+rB) 


(a  2~Hi  i  )  r b 

(al+rB)(a2+rB) 


(40 


P[B] 


(a1-Ka2)(a1-kt2+rB)r2 

ala2(al+rB)(a2+rB) 


r2{  2  (a  i  +a2)  3+l  2  (a  1+^2^  2+a  i +a  ia  2+a2  B+^°  1  +tl2  ^rB^ 

ala2(2al+rB)(2a2+rB)(al+a2+V 


(4  I  ) 


-13- 


'277  -  (U  1.^2  )rB  _  _1 _  8ala2(.aL+a2)'(a?+a2)rB  a2(Sal"rB) 

11 A  '  "  (al+rB)(a2+rB)  (a  2  ~oi  j )  2  (Va2+rB)  (2rtl+rB) 


a  (8a.-r  ) 

+  -k-4  ■?-  ]  +  * 


"V 


(  2a  2  +rB  ) 


.2  2  N  2 

,  (a2"al)rB  "“l11 
{7— t— - -  e 


(.j-Oj)2  a i  f a+rB  ^ 


2  2  2 

(a2“al)rB  ■°2t 
e 


a2(a2+rB} 


a2(3al"rB)rB  _2alt  al(3a2'rB)rB  “2a2t 
+  ttt~ — : - x e  + 7~n ; x - e 


al(2al+rB) 


a2(2a2+rB) 


[6a1a2-(c]-K,2)rB]rB 

~  e  }  . 


(al+“2+rB) 


(4  4) 


ng  (t  )  =  nfi  (t  )'  , 


(45) 


and  again  the  variances  may  be  computed  from  equations  (26)  and  (27). 


5.  COMPARISONS 


The  two  examples  given  above  are  stochastic  formulations  of  the 
Lanchester  Square  Law.  It  will  be  instructive  to  compare  these 
results  with  the  equivalent  Lanchester  results  which  are  given  below 
and  are  very  well  known,  see  for  example,  Clark  [4],  Using  our 
notation, 
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(46) 


nA(t)  **  2  cosh  /rArfi  t 


r'T  slnh  /rArB  c- 
A 


and 


n  (t)  =  cosh  /r  r  t  -  2 J —  sinh  /r.rD  t. 
B  A  B  ✓  r„  A  B 

D 


(47) 


This  formulation  is  in  fact  deterministic  and  we  are  interpreting 
the  variables  as  the  mean  values  of  the  random  variables  NA(t)  and 
Ng(t).  This  is  the  generally  accepted  interpretation  and  indeed  no 
other  viewpoint  seems  reasonable. 

In  making  our  compafisions  we  should  note  that  equation  (46)  has 
the  following  propertied:. 


A  _ 

1.  If  —  <  1/4  , .n  (t)  has  a  zero  at 
rB  A 


and  becomes  negative  thereafter.  Thus  for  t  greater  than  this  value 


we  set  nA(t)  **  0  (we  cannot  have  negative  contestants  on  a  side). 


A  _ 

2.  If  —  >  1/4,  n.(t)  is  always  positive,  but  decreases 
rB  A 

monotonically  from  t  =  0  to 


t 


- -  tanh  ^1/2 


/rArB 
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where  it  has  a  minimum  and  then  starts  to  increase.  At  this  time  we 
continue  the  function  level  at  the  value  it  has  reached  (our  problem 
does  not  allow  for  reinforcements). 


3.  If  — ■  **  1/4,  n  (t)  is  asymptotic  to  zero  as  t  -*• 
rB  A 


n  (t)  has  complementary  properties.  It  goes  to  zero  where 
n  (t)  has  a  minimum;  it  has  a  minimum  where  n  (t)  has  a  zero;  and  it 

A  A 

rA 

is  asymptotic  to  zero  as  t  +  “  for  —  *  1/4. 


The  comparisions  we  are  making  here  are  motivated  by  the  fact 
that  it  is  common  in  combat  models  where  it  is  known  that  the 
interfiring  times  are  not  ned,  to  assume  that  they  are  ned  and  use  the 
means  of  the  true  distributions.  Thus,  if  p  is  the  mean  of  the  true 
interfiring  time  distribution,  the  killing  rate,  r,  is  taken  to  be 
—  (done  appropriately  for  each  side)  and  either  the  "Stochastic 
Manchester"  (both  sides  ned)  or  the  Manchester  equations  are  used  with 

p  ' 

the  appropriate  —  s  as  the  attrition  coefficients. 

The  following  comparative  curves  throw  some  light  on  this  matter, 
at  least  at  the  two-on-one  level.  In  all  of  these  figures  there  are  3 

parameters  p^,  and  r^.  For  comparative  purposes  we  shall  combine 

A 

the  two  paramters  pA  and  p^  wherever  possible  into  —  = 

Pa 

which,  in  fact,  always  occur  together  this  way  in  the  Manchester  (M) 
and  the  "Stochastic  Manchester"  (SL)  (all  ned  stochastic  version) 
equations.  However,  in  our  example  (on  which  the  A  side  has  Erlang 
(2)  firing  times  and  which  we  shall  designate  as  the  E  model)  they  do 


not  in  general,  occur  together.  So  for  comparative  purposes  we  use  pA 
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and  rA  to  characterize  that  situation,  even  though  rA  is  not  the 
instantaneous  kill  rate  (which  is  a  function  of  the  backward 
recurrence  time)  but  simply  the  ratio  of  two  parameters. 

In  Figure  1  we  let  rA  and  rg  be  some  typical  fixed  values,  let  pA 


vary  widely  and  observe  n.(t).  This  makes  no  difference  in  the  L  and 

A 

SL  models  but  causes  considerable  variation  in  the  expected  values 
over  time  in  the  E  model.  We  observe  that  if  rA  is  held  constant  then 
as  pA  increases  and  therefore  pA  increases  (or  — —  the  "firing  rate" 


decreases)  then  the  E  model  n^(t)  decreases.  nA(t)  for  the  E  model 
is  always  less  than  the  SL  model  and  diverges  more  as  pA  increases 
( pA  increases).  The  L  model  may  be  either  greater  than  or  less  than 
either  or  both  the  SL  and  E  models  depending  on  the  values  of  the 
parameters.  The  variations  between  models  may  be  very  substantial 
even  in  the  early  stages  of  the  conflict. 


In  Figure  2  we  have  again  plotted  nA(t)  but  this  time  pA  is 
fixed  and  we  vary  the  parameters  rA  and  rg  such  that  the  ratio 


—  is  fixed.  In  all  cases  n.(®)  is  different  for  each  model  but  is 
rB  A 

the  same  within  each  model  as  the  parameters  vary.  This  is  equivalent 
to  only  varying  the  mean  firing  times  on  each  side,  maintaining  the 
kill  ratio  constant.  What  does  change  is  that  the  time  to  near¬ 
convergence  increases  drastically  as  the  rA's  and  rg’s  decrease.  This 
is  not  unexpected.  And  again  there  is  substantial  variations,  in 


n^(t)  among  the  models. 


In  Figure  3  we  shift  to  plotting  n^(t)  and  hold  the  A  side's 
parameters  constant  and  vary  r^.  In  this  case  the  approach  to  near¬ 


convergence  on  n  ( “)  varies  with  rg.  As  rg  increases  the  approach 


time  decreases,  at  the  same  time  n^i®)  itself  increases.  It  is 
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uniformly  true  that  for  n  (t)  the  E  model  always  is  larger  than  the 

O 

SL  model  and  the  L  model  may  be  located  any  where  relatively  depending 

on  the  value  of  the  parameters.  Again  the  variations  among  the  models 

is  quite  substantial. 

We  note  at  this  point,  that  Clark  [A],  has  previously  shown  that 

the  *L  and  SL  model  vary  and  has  derived  the  error  term  which  he  calls 

•'bias". 

Figure  4  is  similar  to  Figure  1  in  that  rA  and  rg  are  fixed  while 
pA  is  varied,  but  here  we  are  observing  the  standard  deviation  of 
NA(t).  There  is  no  plot  for  the  Lanchester  model  since  it  is  purely 
deterministic.  As  before  varying  pA  does  not  change  anything  in  the 
SL  model  but  the  E  model's  standard  deivation  decreases  as  pA 
increases.  The  difference  may  be  considerable  even  in  the  early 
stages  if  the  kill  rates  are  large.  Here  we  observe  that  the  standard 
deviation  for  the  E  model  may  be  greater  than  or  less  than  the  SL 
model  depending  on  the  kill  rates. 

Figure  5  is  similar  to  Figure  2  as  we  hold  pA  fixed  vary  rA  and 
rB  such  that  rA/rg  remain  fixed  and  observe  the  standard  deviation  of 
NA( t ) .  Again,  as  rA  and  rg  vary  each  model  is  asymptotic  at  infinity 
to  the  same  value  of  the  standard  derivation  but  the  asymptotes  vary 
between  the  SL  and  E  models  although  not  substantially.  The  time  to 
near-convergence  increases  dramatically  as  the  kill  rates  increase. 

Finally,  Figure  6  is  similar  to  Figure  3  where  the  A  side’s 
parameters  are  held  constant,  rg  is  varied  and  the  standared  deviation 
of  Ng(t)  is  observed.  As  before,  the  approach  to  near-convergence  of 
the  standard  deviation  of  Ng(t)  increases  as  rg  increases.  Substantial 
variation  between  models  is  again  evident. 
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6.  Conclusions 


We  have  developed  for  the  first  time  the  general  solution  for  the 
state  probabilities  for  the  two-on— one  stochastic  duel  and  from  these 
derived  the  winning  probabilities  and  the  mean  and  standard  deviation 
functions  for  the  state  random  variables.  These  solutions  are 
illustrated  by  two  examples.  In  the  first,  all  interfiring  times  are 
negative  exponential  (this  is  the  so-called  Stochastic  Lanchester 
version  of  the  two-on-one  Square  Law)  and  in  the  second  example  the 
side  with  two  contestants  had  Erlang  (2)  interfiring  times  each  with 
the  same  means  as  in  example  1. 

These  two  examples  were  then  compared  with  the  corresponding 
Lanchester  Square  Law  solution  to  illustrate  that  all  three  of  these 
models  vary  substantially  from  each  other  and  neither  the  Stochastic 
Lanchester  nor  the  Lanchester  formulation  is  a  satisfactory 
approximation  to  the  general  model  at  least  at  the  two-on-one  level. 
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APPENDIX  I 


In  this  appendix  we  provide  graphs  supplementing  the  ones 
presented  in  the  body  of  the  paper  by  including  more  values  of  the 
parameters  r^,  r^,  and  p^.  The  correspondence  between  the  various 
sets  in  the  text  and  the  appendix  are  shown  in  the  table  below. 


Text 

Appendix 

1 

1 

2 

2.3 

3 

4.5 

4 

6 

5 

7,8 

6 

9  ,  10 

These  figures  demonstrate  again  that  variations  among  the  models 


are  substantial. 
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APPENDIX  II 


For  the  record,  we  give  the  results  of  the  two-on-one  duel  with 
SrLang  (2)  interfiring  time  on  both  sides  here.  Let 


f  X  (XA) 
A 


2  XAe 


T  xBe 


Then  from  reference  [1] 


-  ~  (1  -  /qA  )t  "  IT  (1  +  /qA)C 

fA(0  =  —~r  U  -  e  A  ] 

V,A’^A 


and  f,.(t)  is  the  same  with  all  A's  replaced  by  B's. 


Using  the  notation 
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A  A 
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D  D 
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